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Relativistic bound-state calculations in Light Front Dynamics 
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We calculated bound states in the quantum field theoretical approach. Using the Wick-Cutkosky model and an 
extended version of this model (in which a particle with finite mass is exchanged) we have calculated the bound 
states in the scalar case. 



1. Introduction 

There are several ways of calculating the bound 
states of a system in the context of quantum field 
theory. We have chosen to calculate the bound 
states in a quantum field theoretical approach in 
Light Front Dynamics (LFD), since LFD is very 
well suited for this problem (see e.g. (|]). 

Our goal is to move towards QCD, but before 
treating a realistic system we have done calcula- 
tions in the scalar case. In the case that a mass- 
less particle is exchanged (/i = 0) we treat the 
Wick-Cutkosky model. An extended version of 
this model is used in the case that a particle with 
mass (/i ^ 0) is exchanged. Work on this subject 
has already been done by others; e.g. 

After setting up the equation and explaining 
which method we use to solve it (section 2), we 
will present some results in section 3. At the end 
we will give a short conclusion and outlook. 

2. Setting up the problem and solving it 

The physical system we have looked at consists 
of two scalar particles of equal mass m exchanging 
another scalar particle. This particle can either 
have a mass equal to zero (^ = 0) or it can have a 
mass (i unequal to zero. The Lagrangian for this 
system is given by: 
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where g is the coupling constant. 

Taking along only two and three particle states, 
we can derive the bound-state equation from this 



lagrangian. It is given by: 
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Here M is the total mass of the system and 
K{p^,x]p'^ 1 x') is the kernel of the equation. 
There are several possible expressions for the 
kernel, depending on what approximation we 
make. Here we have been working in the lad- 
der approximation. In first order we have two 
time ordered diagrams (Fig. |l|). Each of these 
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Figure 1. First order time-ordered diagrams in 
the ladder approximation for one particle ex- 
change. 

diagrams corresponds to an energy denominator 
which is present in the kernel of the bound-state 
equation (Eq. (||)). The total kernel contains, 
besides the energy denominators, a phase-space 
factor and a theta function, which determines in 
which time-ordering we are working. 
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For the total kernel we can write: 
K(p ± ,x; p' ± ,x r ) = 

1 / 9(x'—x) . 6(x' — x) 
<S>(x,x') \2(x> -x)D a + 2(x-x')D b 



(3) 



Here $(a;, x 1 ) is representing the phase-space fac- 
tor, which is given by: 



$(x,x') = y/ x{l - x)x' {1 - X') , 



(4) 



and D a and Db are the energy denominators cor- 
responding to left or right graph in Fig. |l|. The 
expressions for these energy denominators are: 
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The bound-state equation has been solved by 
making an expansion into basis functions, sym- 
metrizing it and then using Gaussian integration 
to solve it. The wave functions depend on both x 
and \p\ and we have used two different functions 
to expand in. For the x part of the wave function 
we have used cubic spline functions and for the 
\p\ part we have used a basis which was first in- 
troduced by Olsson |(| and Weniger M . We have 
tested the latter basis in instant form H and we 
have seen that it works well. 

3. Results 

Using the method described in the previous sec- 
tion, we are able to calculate the spectra and wave 
functions in case of massless exchange as well as 
in the case of massive exchange. Here we have 
to remark that we did not take into account any 
higher order terms in the Hamiltonian or self en- 
ergy diagrams. 

First of all we have looked at how many func- 
tions we need to take along to reach a reason- 
able convergence. A reasonable number of spline 
functions to take into account in the calculations 
is 10. Although 8 spline functions is a bit on the 



low side, it already gives good results as well. For 
the Olsson Weniger functions we need to take 3 
functions to reach a good convergence. Taking 
less functions already gives a good convergence 
in the lower states, but not in the higher excited 
states. Note that increasing the number of Ols- 
son Weniger functions in the expansion causes a 
large increase in computer time. Due to this we 
have chosen to use 8 spline functions and 3 Ols- 
son Weniger functions in the expansion to obtain 
the results presented here. 

In the case of massless exchange (fi — 0.0) we 
have calculated the masses given in the second 
column of table [l]. In the third column of this 
table the masses of the bound states are given 
in the case of massive exchange (/i = 0.15). In 
both cases we have used a mass m = 1.0 for the 
particles and a coupling constant of g = 17.36. 



n 


mass (/i = 0.0) 


mass (/i = 0.15) 





0.75 


0.91 


1 


1.63 


1.75 


2 


1.74 


1.86 


3 


1.86 




4 


1.91 





Table 1 

Masses corresponding to the bound states cal- 
culated in the case of massless and massive ex- 
change using the following parameters: m = 1.0 
and g — 17.36. 

We have also calculated the wave functions cor- 
responding to the masses given in table [l]. For 
the massless case the wave functions correspond- 
ing to the first three states are given in Fig. |^, 
|^ and The wave function corresponding to 
the ground state in the case of massive exchange 
(/i = 0.15) is given in Fig. || Looking at these 
figures one should notice that the wave functions 
are not drawn in the ordinary way. This is due to 
the fact that the wave function depends on x (the 
fraction of the p + -momentum) and pj_, instead of 
p as is usually the case. 

4. Conclusions and outlook 

The bound-state equation in the scalar case is 
solved by making an expansion into basis func- 





Figure 2. Ground state wave function for massless 
exchange. 




Figure 5. Ground state wave function in the case 
of massive exchange (fi = 0.15). 



tions. The basis we have used is well suited for 
this problem and we expect it to work in more 
realistic cases as well. With this method we can 
calculate the ground state as well as excited states 
in both cases, massless and massive exchange. 

The plans for the future are to move towards a 
more realistic system, i.e. to include spin. And 
also to included higher order terms in the Hamil- 
tonian and see whether these terms give a large 
contribution. 



Figure 3. Wave function of the first excited state 
in the case of massless exchange. 
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Figure 4. Wave function of the second excited 
state in the case of massless exchange. 



